Astronomy <fe Astrophysics manuscript no. paper'dyna'2 


© ESO 2008 


June 26, 2008 





Gas flow models in the Milky Way embedded bars 

N. J. Rodriguez-Fernandez^ and F. Combes^ 

1 IRAM, 300 me de la Piscine, 38406 St. Martin d'Heres, France 

e-mail: rodriguez@iram.fr 
^ LERMA, Observatoire de Paris, 61 Av de I'Observatoire, 75014 Paris, France 

Received September 15, 1996; accepted March 16, 1997 

ABSTRACT 

Context. The gas distribution and dynamics in the inner Galaxy present many unknowns as the origin of the asymmetry of the 
Zw-diagram of the Central Molecular Zone (CMZ). On the other hand, there are recent evidences in the stellar component of 
the presence of a nuclear bar that could be slightly lopsided. 

Aims. Our goal is to characterize the nuclear bar observed in 2MASS maps and to study the gas dynamics in the inner Milky 
Way taking into account this secondary bar. 

Methods. We have derived a realistic mass distribution by fitting the 2MASS star counts map of Alard (2001) with a model 
including three components (disk, bulge and nuclear bar) and we have simulated the gas dynamics, in the deduced gravitational 
potential, using a sticky-particles code. 

Results. Our simulations of the gas dynamics reproduce successfully the main characteristics of the Milky Way for a bulge 
orientation of 20° — 35° with respect to the Sun-Galactic Center (GC) line and a pattern speed of 30-40 km s~^ kpc~^. In our 
models the Galactic Molecular Ring (GMR) is not an actual ring but the inner parts of the spiral arms, while the 3-kpc arm and 
its far side counterpart are lateral arms that contour the bar. Our simulations reproduce, for the first time, the parallelogram 
shape of the to-diagram of the CMZ as the gas response to the nuclear bar. This bar should be oriented by an angle of 
~ 60° — 75° with respect to the Sun-GC line and its mass amounts to (2 — 5.5) 10^ Mq. We show that the observed asymmetry 
of the CMZ cannot be due to lopsidedness of the nuclear bar as suggested by the 2MASS maps. 

Conclusions. We do not find clear evidences of lopsidedness in the stellar potential. We propose that the observed asymmetry of 
the central gas layer can be due to the infalling of gas into the CMZ in the 1=1.3° -complex. 

Key words. Galaxy: structure - Galaxy: center - Galaxy: kinematics and dynamics - ISM: kinematics and dynamics - Methods: 
numerical 



1. Introduction 

Studying the structure and the dynamics of the Milky Way 
galaxy is necessary to understand its formation and evolu- 
tion. For example, the large scale star formation rate, the 
gas transfer to the nucleus and the fueling of the central 
super-massive black hole are strongly influenced by the 
large scale dynamics of the Galaxy. 

At least three components are needed to explain the 
rotation curve of the Galaxy: a disc, a bulge and a mas- 
sive dark halo. In near-infrared images, the bulge is a tri- 
axial structure difficult to distinguish from a bar. The 
first direct evidence of the barred nature of the Milky 



Way was found by Blitz & Spergel (1991). The analysis 



of COBE/DIRBE data ( Dwek et al 



1997 



suits. 



Dwek et al. 



1995 Binney et al. 



Freudenreichj 1998') supported Blitz & Spergel re- 
(1995) found that the peanut-shaped 



bulge seen in the COBE image is best fitted by boxy 
Send offprint requests to: N. J. Rodriguez-Fernandez 



Gaussians functions. Since COBE, the direct evidence of 
the bar has increased based on photometric surveys and 



star counts in the near-infrared (see for instance Lopez- 
ICorredoira et al.||2000| |2005| [Babusiaux fc Gilmore||2005| 
and references therein). Although, there is a consensus 
on the triaxial nature of the bulge, the actual shape 
and the inclination is still under debate with currently 
no clear agreement among different authors. The inclina- 
tion angles found range from ~ 10° — 15° (Binney et al. 



1991 Freudenreich 1998 Lopez-Corredoira et al. 1997) 



to ~ 30° - 45° (Weiner & Sellwood 1999 Sevenster 1999 



Lopez-Corredoira et al.|2001 ) with a number of works giv 
ing values in the range 15° — 30° (Mulder & Liem [1986 



Binney et al.|1997| IFux"1999VEnglmaier fc Gerhard'1999 
Bissantz fc Gerhard 2002; Bissantz et al. 2003: Babusiaux 
fc Gilmore||2005| |L6pez-Corredoira et al.||2005| . 



In addition to the bar-like boxy-hulge, some recent 
works have also proposed the existence of another bar that 



would be longer and thinner than the bulge ( Picaud et al 
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2003 



|L6pez-Corredoira et al.|[200T| [20071 Pelij amin et al.| To our knowledge, the wor ks by|Binney et al.| ( |1991[ ) 



2005 1. As the bulge, the near side of this long and thin bar 
will also be located in the first Galactic quadrant. In con- 
trast, the long bar would be oriented at an angle of ~ 45° 
with respect to the line of sight towards Galactic center. 
In the present paper we will not discuss this component 
as the mass should be dominated by the bulge. Hereafter, 
we will refer to the triaxial boxy bulge as "large bar" or 
simply "bulge" . 

Due to the difficulty to measure the stellar veloci- 
ties, most of our knowledge of the Galactic dynamics 
come from spectroscopic observations of the atomic and 



Jenkins & Binney ( 1994 1 and Lee et al. (|1999| are the 
only dynamical models that have tried to interpret the 
dynamics of the complex Galactic center region (the cen- 
tral 500 pc) . Other models of the Galactic dynamics have 
been mainly devoted to larger scales, from the works of 



Mulder fc Liem| (|1986|) and Wada et al (1994) to the more 
recent works by Fux (1999), Englmaier & Gerhard (1999|, 



Weiner fc Sellwood] ( |1999f and |Bissantz et al.| ( |2003| ). The 
bar pattern rotation speed for most of those models are 
in the range from 40 to 60 km s~^ kpc~^. Regarding the 
corotation radius, it is estimated to be in the range 3-4 



molecular interstellar medium (Rougoor & Oort 


1960 


Scoville 


19721 


Liszt & BurtonI |1978| Bally et al. 


1987 


Dame et al.||2001 


). The gas exhibit large non-circular mo- 



kpc by some works ( Englmaier fc Gerhard|[T999 Bissantz 
et al.|[2003 Habing et al. 20061 and in the range 4-5 kpc 



tions in the inner Galaxy. Among the possible explana- 
tions for these kinematics, de Vaucouleurs (1964) proposed 
the gas response to a bar potential. This explanation has 
also been discussed in other observational works as those 



by other wo rks (|Combes|1996[ |Sevenster|1999{ |Fux|1999 1 

In contrast, Binney et al. (1991) placed the corotation at 
2.4 kpc and found a pattern speed of 80 km s~^ kpc~^. 
As pointed out by ICombesI (|1996|) , the discrepancy be- 



of Kerr (1967), Bania (1977) or Liszt fc Burton (1978) 



Nevertheless, it is the comparison of the observations with 
numerical simulations of the gas dynamics what have al- 
lowed to better understand the structure of the Galaxy 
and to derive important dynamical information as the pat- 
tern speed of the bar. The first attempt was that of [Mulder] 
fc Liem ( 1986 1 solving the gas-dynamical equations to get 



tween Binney et al. estimation and other estimations can 
be due to the fact that the molecular gas distribution in 
the Galactic center region used by Binney et al. could be 
the response to a small nuclear bar instead of the large 
bar of the Galaxy. 

Indeed, there are recent evidences of the presence of 
a distinct structure in the nuclear region of the Galaxy. 



the gas flow in a weak bar for the whole Galaxy. 

A major step in our understanding of the gas kinemat- 
ics in the Galactic Center was the paper by [Binney et al.| 
(1991 1. These authors made a comparison of the CO(l-O) 
and HI longitude- velocity diagrams (hereafter Zu-diagram) 
of [Bally et al.|([1987[) andjBurton fc Liszt[([1978[), with the 



The star counts map made by Alard (2001) with 2MASS 
data shows an excess of counts in the inner four degrees of 
the Galaxy that resembles a small nuclear bar. This bar 
could be lopsided with the barycenter displaced to nega- 



tive longitudes. Nishiyama et al. (2005) have also found 



Iv locus of closed stellar orbits in a bar potential (xl and 
x2 orbits). The Zti-diagram of the molecular gas in the in- 
ner ~ 4° of the Galaxy shows a parallelogram shape that 
resembles the Iv locus of the cusped xl orbit for a bar with 



evidence of a structure different to the large bar or bulge 
in the inner degrees of the Galaxy. 

In the present paper, we want to investigate the gas re- 
sponse to the mass distribution inferred from the 2MASS 
data. In particular, we model the central mass distribution 
as a small nuclear bar. We are interested in the role of this 
nuclear bar to explain the differences among the results of 



metry of the CO parallelogram, which is lopsided towards 
positive longitudes while the modeled parallelograms are 
centered at longitude 0°. In addition the expected upper 
right vertex of the locus of the cusped xl orbit has a veloc- 
ity close to 0, while the upper right vertex of the observed 



inclination of ~ 16°. In spite of the globally successful in 
terpretation of the CO and HI Zu-diagrams, the [Binney[ previous works (Binney et al. 1991 Fux 1999 Combes 
et al.[([1991[) model cannot account for the observed asym- 1996). We also want to study whether a lopsided nuclear 

bar could explain the apparent lopsidedness of the molec- 
ular gas distribution in the nuclear region, whose com- 
plex dynamics are far from being understood. Therefore, 
we have modeled the 2MASS star counts maps of Alard 
fitting different parametric functions using three compo- 
nents: disc, bulge/bar and nuclear bar. We use that mass 
distribution to determine the galactic stellar potential and 
to model the gas dynamics using a N-body code that takes 
into account cloud collisions. 

Nested bars are common in external galaxies. As much 
as one third of barred galaxies have a secondary nuclear 



CO parallelogram shown in Fig. 2 of Binney et al. ( 1991 1 
has velocity of ~ 150 km s^^. 

The qualitative insight on the gas dynamics that one 
could obtain by comparing the Zti-diagram of the gas with 
that of a few stellar orbits, should be tested performing 
simulations of the gas dynamics that take into account 
the dissipative nature of the gas clouds. This was done 



by Jenkins fc Binney (1994), who tried to reproduce the 



/?;-diagram of the CMZ using a sticky-particles code but 
with moderate success. Using a combination of SPH and 



bar (Laine et al. 2002). In some cases the secondary bar 
is coupled with the primary bar ( Shaw et al.[[1993 ) while 
in other cases the two bars rotate with different speeds 
(Wozniak et al. 1995). The present paper is the first one 



coUisionless particles, Lee et al. (1999) were able to re- 



produce an inner ring with the right inclination in the Iv- 
diagram and also some of the gas clouds with non-circular 
velocities in the range / < |6°|. 



dealing with the gas flow in the Milky Way nested bars 
and we will assume that both bars are dynamically cou- 
pled (they rotate with the same speed). As we discuss in 
Sect. 8, this seems a good assumption since our Simula- 
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tions explain for the first time some characteristics of the 
inner Galaxy. 

The paper is organized as follows: in the next section 
we summarize our understanding of the distribution and 
kinematics of the interstellar gas. In Sects. 3 and 4 we 
present the 2MASS star counts model and the fits to the 
data. The gas dynamics simulations are described in Sect. 
5, and the results are presented in Sects. 6 an 7. We discuss 
the nature of the nuclear bar in Sect. 8. Our conclusions 
are summarized in Sect. 9. 

2. The distribution and kinematics of the 
interstellar gas 

In this section we summarize outstanding structures in the 
longitude-velocity diagram (hereafter /u-diagram) of the 
interstellar gas. The large scale ??;-diagram derived from 
observations of the CO(l-O) line by ( Dame et al.||2001 l is 
shown in Fig. [l] while Fig. |2] shows the /z;-diagram of the 
nuclear region derived also from CO(l-O) observations by 
Bally et aT] ( |1987[ ) . We also discuss the plausible face on 
view of the Galaxy. The different features discussed here 
will be used in the following sections to compare with our 
simulations of the gas dynamics. 



2.1. The terminal velocity curve 

The terminal velocity curve (TVC) is not a kinematic fea- 
ture in itself as many independent features contribute to 
this curve at different Galactic longitudes. However, we 
discuss it here since it has been used to constrain several 



dynamical models of the Milky Way ( Binney et al. 


1991 


Weiner k. Sellwood 


|1999| lEnglmaier & Gerhard 


1999 


Bissantz et al. 2003 


1. There has been different determi- 


nations of the TVC using CO and HI data (Fich et al. 


1989 


Burton & Liszt 1993| Clemens 1985| Alvarez et al. 


1990 


1. The dispersion on the TVC points derived by dif- 



ferent authors is typically 10 km s (see for instance 



Fig. 6 of Bissantz et al. 2003 ) . The TVC measurements by 
Fich et al.| ( |1989[ ) and Burton fc Liszt] ( |1993[ ) are shown 
in Fig. [1] and in the figures showing the simulations re- 
sults. In addition, in Fig. [l] we have indicated with boxes 
the Sagittarius, Scutum, Norma and Centaurus tangent 
points, located respectively at I ~ 50°, 30°, —30° and 
—50° (see for instance Dame et al 
Gerhard|[l999| . 



2001 Englmaier fc 



2.2. Spiral arms 

Many works have studied the spiral structure of the 
Galaxy using different tracers as HII regions, atomic and 
molecular gas or photometric observations in different 
bands from the near to the far infrared. The different mod- 



els of the spiral arms have two (|Bash ( 1981 1 using HII 
Drimmel fc Spergel ( |2001[ ) using K-band observa- 



regions, 
tions. 



Benjamin (20081 using Spitzer), three (Nakanishi fc 



Sofue 



2003 Kulkarni et al. 1982 using HI) or four arms 



fc Sofue] p006| using CO or ]Drimmel] ( ]2000[ ) using the 240 
microns luminosity). Deprojecting a given data set to find 
out the face on view of the Galaxy is a difficult task in 
particular for the inner Galaxy. For instance, the kinemat- 
ical distances based in the rotation curve are not accurate 
due to the presence of important non-circular motions. 
Furthermore, after performing the deprojection there is 
an uncertainty in how different structures are associated 
to a common arm (see [Nakanishi fc Sofue|2006] ). Different 
authors have even fitted different spirals to the same data. 
A plausible explanation for the discrepancy among the two 
arms derived from stellar tracers (K band) and the four 
arms that are found in most of the interstellar tracers is 
that the Milky Way has two stellar spiral arms but the gas 
response is a four-armed spiral (Drimmel 20001. In any 



case, the major arms are logarithmic with pitch angles of 
10° — 15°. In addition, to these grand-design spirals there 
is a more flocculent structure giving a number of minor 
arms like the Local Arm, which is a small arm linked to 
the Perseus spiral arm. In order to compare with our sim- 
ulations, in Figs. |8l ]9land]T0l we have traced the sketch of 



the logarithmic spiral arms of]Nakanishi fc Sofue] (]2006|, 
which are very similar to the classical spirals of Georgelin] 
fc Georgelin (19761. The locus of different spiral arms in 



the Zu-diagram are shown in Fig. ]T] The Carina arms de- 
serves a special mention since the Carina tangent point is 
clearly seen at I ~ —80°. 

The determinations of the speed of the spiral pattern 
give values in the range from 13.5 to 59 km s~^ kpc~^ 



(Lin et al. 


1969 Fernandez et al. 


2001 Bissantz et al. 


2003 


Martos et al. 


2004 


Debattista et al. 


2002 


). Thus, 



the spiral pattern speed could be lower that the bar one. 



2.3. The Galactic Molecular ring 

The so-called Galactic Molecular Ring (GMR) is one of 
the most prominent structures in the Iv-diagram (Fig. [T]) . 
The locus of the GMR is more or less a diagonal line that 
passes trough the origin as a circular ring would do. If the 
GMR is a real ring, it could be associated to a resonance. 



Binney et al. ( 1991 1 proposed that the GMR could be 
located at the Outer Linblad Resonance. Other models 
that place the corotation farther out have suggested that 
the GMR could be better explained by the Ultra Harmonic 
Resonance ( Combes]] 1996 ). 

However, the real spatial distribution of the GMR is 
not known and it could be composed by imbricated spiral 
arms instead of being an actual ring. Some authors (see 
for instance Englmaier fc Gerhard 19991 consider that. 



if formed by several arms, the GMR should be located 
outside of the corotation since they exhibit almost circu- 
lar velocities. In their face-on deprojection of the neutral 



and molecular gas distribution, Nakanishi fc Sofue (20061 



(Georgelin fc Georgelin ( 1976[ ) using HI regions; Nakanishi 



found that the GMR seems to be the inner part of the Sgr- 
Carina arm (the far-side of the GMR) and the Scutum- 
Crux arm (the near-side of the GMR) . 
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Fig. 1. Longitude-velocity (Iv) diagram of the CO(l-O) emission ( Dame et al.|2001 l. The solid lines trace the position 
of some remarkable features as the locus of the spiral arms, the 3-kpc arm and the Connecting Arm. The black dashed 
lines indicate the contour of the Galactic Molecular Ring. The solid circles are the terminal velocities measurements 



of Fich et al. ( 1989 1 using CO while the triangles are the terminal velocities determined from the HI data of Burton & 



Liszt (19931. The boxes mark the position of the Sagittarius, Scutum, Norma and Centaurus tangent points, located 
respectively at I ~ 50°, 30°, —30° and —50°. The lines concerning the Nuclear Disk or Central Molecular Zone are 
blue but they are better shown in the next figure. All these lines can be used to compare with the figures showing 
the simulations results. 



2.4. Inner arms: the 3-kpc arm 

The 3-kpc arm is clearly seen in the Zw-diagram of Fig. 
[l] with a velocity of -53 kms'^ at I = 0°. The name of 
this feature comes from a tangent point at I — —22° that 
corresponds to a distance of ~ 3 kpc. This arm is located 
in between the Sun and the GC since the clouds in the 
arm absorbs the radiation from the continuum sources in 
the GC. The 3-kpc arm cannot be explained with a com- 
bination of rotation and expansion ( [Burke fc Til^[l964| 
and has been interpreted as a stationary density wave in a 
barred potential ( Mulder fc Liem||1986 l. The simulations 
of Fux (19991 suggest that it is a lateral arm that sur- 



rounds the bar while, in the simulations by Englmaier fc 



Gerhard ( 1999 1 , it would be an small arm arising from the 



extremity of the bar. 



Sevenster (19991 and Habing et al. (2006) propose a 



different explanation. Sevenster ( 1999 1 found old stars as 



sociated with the 3-kpc arm. They interpret the fact that 
both old stars and gas follow the same trajectories as the 
probe that the arm would be the locus of closed orbits and 



not a spiral density wave maximum. They propose that the 
3-kpc arm has its origin near one of the two points where 
the bar meets its corotation radius and that the arm can 
be a channel to transport gas from the corotation to the 
GC and to fuel the star formation in cloud complexes like 
Sgr B2. 

The recent discovery of the far-side counterpart of the 



3-kpc arm by Dame & Thaddeus (2008 1 will certainly con- 



tribute to understand the exact nature of these features 
(see Sect. 6.3). 

2.5. The connecting arm and Bania's clumps : tracing 
the dust lanes ? 

The "Connecting Arm" is clearly seen in Fig. [T] as a fea- 
ture with high inclination with a velocity of ~ 100 kms~^ 
at I = 10° and more than 250 kms~^ at Z = 3° . The name 
"Connecting Arm" comes from the fact that it seems to be 
connected to the nuclear ring. This a very puzzling struc- 
ture that has not been discussed much in the literature 



(see Fux 1999 for a short summary). On the one hand. 
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Burton & Liszt ( 1978 1 and Liszt & Burton ( 1980 1 sug- 



gested that the Connecting arm could be the edge of the 
HI ring found at a radius of ~ 1 kpc (see below). On the 



zoo 



other hand, in one of his simulations, Fux| ( 1999 ) obtained 
a feature in the /t^-diagram that resembles the Connecting 
Arm. In that simulation, the Connecting Arm-like feature 
is the locus of the off-axis shocks or dustlanes (since in 
the optical images of barred galaxies the off-axis shocks 



are seen as dust lanes) . Marshall et al. ( 2008 1 also inter- 
pret the Connecting Arm as the near side dustlane. 

However, in other of the simulations discussed by |Fux 
( |199 9l) the locus of the dustlane in the /u-diagram is a ver- 
tical feature extending approximately from a velocity of 
kms^^ to 200 kms^^ at an almost constant longitude 
(see his figures 15 and 16). This reminds the "clumps" 
discussed by Bania ( 1977[ ). These clumps are outstand- 
ing cloud complexes found at I — 5.5° and I = 3.1° that 
exhibit a huge velocity dispersion of ^ 200 kms^^ in a 
very narrow range of longitude. A portion of the clump 
located at / = 3.1° (known as Clump 2) is showed in Fig. 
2) The Clump 2 has been interpreted as a dustlane by 
Stark & Bania ( 1986) , while in the Fux ( 1999 1 interpreta- 
tion, Banias's clumps would be gas clouds that are about 
to enter the dustlanes shocks, although the actual dust- 
lane would be the Connecting Arm as mentioned above. 



This is also the interpretation of Liszt (2006 2008), who 
have studied in detail the spatial and velocity structure of 
all these features. 



2.6. The HI ring 

The HI emission in the GC region can be modeled as a 
tilted circular ring with inner radius of ~ 300 pc and outer 
radius of 1-1.5 kpc ( Burton fc Liszt|1978 1. This first model 
was improved by Liszt & Burton (1980), who discussed 



that an elongated ring could also explain the observations. 
The spatial distribution of the interstellar gas in the in- 
nermost 3 kpc of the Galaxy has recently been reviewed 
by Ferriere et al. (2007). The model that they consider 
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Fig. 2. Upper panel: Longitude-velocity (Iv) diagram of 
the CO(l-O) emission in the Central Molecular Zone us- 
ing data of Bally et al. ( 1987 ) and integrating all negative 



Galactic latitudes. The solid lines trace the position of 
some remarkable features as the Clump 2, or structures K 



and J of Rodriguez- Fernandez et al. (2006). The dashed 



line indicates the contour of the very crowded central re- 
gion. Lower panel: Integrated intensity map. 



counterpart. Figure |2] also shows the /w-diagram of the 
CMZ. This plot has been obtained integrating the C0(1- 
0) data of Bally et aT] ( |1987i ) for negative Galactic longi- 
tudes to remark a number of kinematical component as 
the Clump 2 (at Z = 3.1°), and the features K and J dis- 



cussed by Rodriguez-Fernandez et al. ( 2006 1 



The kinematics of the inner region are very complex 
with a high fraction of the gas exhibiting non-circular ve- 
locities. Figure [2] shows with a dashed line the contour of 
the inner CMZ that contains the feature M of [Rodriguez- 1 
Fernandez et al.] ( [2006 1 and the inner arms found by 
Sofue (19951. This contour resembles a parallelogram but 



more plausible, despite its inherent uncertainties, is that (figgij) work 



it is not the parallelogram that inspired the Binney et al 



of Liszt & Burton ( 1980 1 since an elongated ring seems 
to be easier to understand if it is composed by gas clouds 
moving in elongated xl orbits. Therefore, they modeled 
the HI ring as an elliptical ring with semi-major axis of 
1.6 kpc, and axis ratio of 3.1, with an inner hole of semi- 
axis 800 pc X 258 pc that encloses the Central Molecular 
Zone (see bellow). 

2. 7. The Central Molecular Zone 

The Central Molecular Zone (hereafter CMZ) refers to 
the central accumulation of gas in the inner hundreds of 
parsec of the Galaxy. This gas is mainly molecular and ex- 
tends continuously approximately from —1.5° to 2° (see 
Fig. |2| . The spatial distribution is not symmetric due to 
the prominent cloud complex located at / = 1.3° (here- 
after 1=1.3° -complex) that lacks of a negative longitude 



The whole Zw-diagram shown in Fig. [2| is very dif- 



ferent to that used by Binney et al. (1991). These au- 



thors integrated the CO data of Bally et al. ( 1987 1 over 
positive and negative latitudes and studied the region 
2.2° > I > -2.2°. In their Fig. 2, the feature K is not 
clearly detached from the inner Iv parallelogram and they 
considered that it is part of a larger Iv parallelogram. 



However, following the analysis of the Bally et al. ( 1987 1 



data by Rodriguez-Fernandez et al.l ( 2006 1 , we reckon that 



feature K is a distinct kinematic feature. This conclusion is 
also supported by the large scale (13° > I > —12°) C0(1- 
0) survey by Bitran et al. (1997), whose Zw-diagram shows 
clearly that the feature K extends to Z ~ 3° where it is 
connected to the Connecting Arm. 

Interestingly, when one compares the Iv parallelograms 



of the xl orbits displayed in Fig. 3 of Binney et al. ( 1991 1 
to the parallelogram shown in our Fig. |2| the upper right 
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vertex are in good agreement, showing velocities of ^ 
km s~^. Thus, one of the apparent discrepancies (see Sect. 



1) between the Binney et al. ( 1991 1 model and the CO Iv- 



diagram disappears. 



3. Modeling the star counts map 

3.1. The star counts map 



We have modeled the 2MASS star counts map of |Alard| 
( 2001 1 to determine a realistic stellar potential. In this sec- 



tion we summarize the methods and the main conclusions 
of Alard's work. 

Using the H and K bands and assuming that we know 
the intrinsic color of the sources {H and K), one can cor- 
rect for extinction one of the bands as follows: 

m^^ ^TUK^ [Ak/{Ah - Ak)] [{mn -tuk)- {H - K)]{\) 

where ttik.h are the observed apparent magnitudes in 
K and H band (thus suffering extinction), Ak,h is the 
extinction in each band and m'^^^ is the extinction cor- 
rected magnitude in K band. At the distance of the bulge 
a cutoff rriK = 9 (see below) implies that the sources 
are early M-giants (see Lopez-Corredoira et al. 2001). 
Assuming an intrinsic color iJ — if for M giants of 0.17 
( Binney fc Merrifield|1998||Wainscoat et al.|199"2| and tak- 
ing into account the relative extinctions Ak/Ay = 0.112 
and Ah/Ay = 0.175 ( |Rieke fc Lebofsky|[l985| , one can 
write: 



TO 



K 



rriK — 1.77 {rriH — rnx) + 0.30 



(2) 



Alard ( 2001 1 defined an extinction corrected magni- 
0.30 as 



^corr 
''K 



tude TOe = TO 

TOg = TO/f — 1.77 {rriH — to^) 



(3) 



and constructed his star counts map using this magni- 
tude with a cutoff of 9 mag. In most of this map, the 
density profile at constant longitude is exponential. This 
exponential profile is also present in numerical simula- 
tions. Combes et al. ( 1990 1 showed that a disk with a 



small bulge in its center can evolve in a peanut shaped 
bar with a nearly exponential profile perpendicular to 



the plane of the galaxy. However, Alard (2001) noticed 



that the observed density profile in the inner four degrees 
of the Galaxy shows an excess of sources in the plane. 
Subtracting an exponential profile corresponding to disk 
and bulge, he found that the residual seems to be a small, 
lopsided, nuclear bar in the galactic center. This struc- 
ture has not been found in previous studies since the data 
lacked either the depth or the resolution. 

3.2. The model 

In this section we describe how we have modeled the 
2MASS star counts map. The number of stars in a vol- 
ume dV located at x with a magnitude in the interval M 
and M + dM is: 



dN = $(Af, x)dMdV 



(4) 



If the distribution of magnitudes is the same in all 
the space we can write <i> as <i>(M, x) = (j>{M)p{x), where 
(j){M) is the luminosity function (the fraction of stars with 
magnitude M) and p(x) is the star density in the point 
X. Therefore, the total number of stars with a magnitude 
M lower than a cutoff Air is: 



N = 



dN = 



(l){M)dMp{x,y,z)dV (5) 



As K band luminosity function (j){M), we have used 
the function given by Wainscoat et al.|( 1992l, which is 
very similar to that of Eaton et al.| (1984). However, we 
work with extinction-corrected apparent magnitudes to, 
which are related to the intrinsic magnitudes by: 



M : 



5-5 logio 



(6) 



where s is the distance in pc. Therefore we actually deal 
with a function of the apparent magnitude and the dis- 
tance 4>'{rn, s). 

Our model of the density distribution is given by three 
components: a triaxial bulge or large bar, an exponential 
disk and a small nuclear bar. The exact functional form 
of each component is given below. 

3.2.1. Disc 



We use the exponential disc of Wainscoat et al. ( 1992 1 
which is defined as: 



p (r, z) = po exp - 



{t-DI_\z\ 

,j. lip. 



(7) 



where po stars number density in the solar neighbor- 
hood, D is the distance from the sun to Galactic center, 
hz and hj. are a vertical and radial scale parameters, re- 
spectively. 

3.2.2. Triaxial bulge 

From the COBE/DIRBE images we know the bulge of 
the Milky Way is boxy ( |Dwek et al.][l995|. The 2MAS S 
data confirm this result (Lopez-Corredoira et al. 2005). 



Therefore, we have adopted a hoxy Gaussian (Dwek et al. 
function G2) to represent the bulge. 



p^{x ,y ,z) = poexp(-0.5r^) 

2 



1/4 



(8) 
(9) 



The bulge function is expressed in a coordinate sys- 
tem (x',y',z') that follows its symmetry axes. Therefore, 
a rotation is needed to get the expression in the (x,y,z) 
coordinate system. 



X cos + 
~x sin 
z 



y sin f3^ 

f y cos (3^ 



(10) 



N. J. Rodriguez-Fernandez and F. Combes: Gas flow models in the Milky Way embedded bars 



7 




Fig. 3. Sketch of the different coordinates systems and 
angles used to define the bars in the Galactic plane. The 
z and z' axis are coincident and the 2 > semiaxis is in 
the northern Galactic hemisphere. 



Figure [3] shows the (a:, y) and (a:', j/') axis seen from the 
z = z' > hemisphere. The figure also shows the defini- 
tion of the angle /3^, which is measured counterclockwise 
from the x > semi-axis. However, many papers measure 
the inclination of the bar clockwise from the Sun-GC line 
to the near side of the bar. Therefore, for an easy compar- 
ison with previous results we have defined the angle oP 
as shown in the figure. 



3.2.3. Triaxial nuclear bar 

In contrast to the bulge, the small nuclear bar found by 
Alard does not seem to be boxy. Therefore we have taken 
a triaxial Gaussian function (Dwek et al. function Gl) to 
represent this component: 



p^(a;",2/",2") = Poexp(-0.5r^) 



1/2 



(11) 
(12) 



Where x\ allows to fit a lopsided bar. To obtain the 
bar density in the {x,y,z) coordinate system one should 
apply the rotation: 



3.2.4. Final star counts equation 

Finally, we must also convert p(x,y,z) to p'{l,b,s) since 
we want to model the star counts in Galactic coordinates 
(I, b) and to integrate along the line of sight (s). If {x, y, z) 
is a right-handed reference frame centered in the galactic 
center, with the plane of the sky in the a;z-plane with a; > 
for I < 0, and 2 > towards the Galactic north pole, and 
assuming that the sun is located in {x,y,z) = {0,—D,0) 
then: 



X — —s cos(fe) sin(/) 

y = —D + s cos{b)cos{l) 

z = s sin(6) 



(14) 



The star counts equation towards (Zqi^o) is then: 
N{lo,bo,mc) = AlAbx 

X /^'""^ dm /q'™""' ds (j)'{m, s) p'{Iq, &07 *) cos b 

where cos b is the Jacobian of the coordinates trans 



(15) 



formation (Eq. 14 1. We have assumed that the density is 
constant in the interval from Iq to + Al and from bo to 
bo + Ab. 

We have computed the integral in magnitudes from 
^min — —10.0 mag to rrijnax — 9.0 mag with Am — 
0.5 mag. The integral on the distance along the line of 
sight s has been computed up to Smax — 20 kpc using a 
variable step As: 

0.002 kpc s,R< 0.100 kpc 

0.008 kpc s,R< 0.300 kpc 

0.010 kpc s,i? < 0.700 kpc 

As = ( 0.050 kpc s, i? < 1.500 kpc (16) 

0.100 kpc s,R< 3.100 kpc 

0.200 kpc s,i? < 6.300 kpc 

0.300 kpc s,i? > 6.300 kpc 



4. Fitting the star counts map 

We have fitted the star counts map using the model 
described above. This model has 14 free parame- 
ters corresponding to the disc (po, h^, h^), the bulge 
iPo,xo,yo,zo,a) and the nuclear bar {po,xo,yo,zo,a,xi). 

The star counts equation is evaluated for an initial 
guess of the density function parameters. A parameter 
is computed as the difference of the star counts data points 
{NdataJ and the model (Nmodeh)- 



x" — X cos P'' + y sin 
y" — —X sin fi^ + y cos 
z" = z 



(13) 



The definition of (3^ is exactly analogous to that of 
shown in Fig. [3] Again, in the following, and for coherence 
with most of the literature in the subject, we will use an 
angle measured clockwise from the Sun-GC line to refer 
to the inclination of the nuclear bar. 



1 



data 



N„ 



odel 



i=l 



(17) 



where n is the number of points to be fitted and a is 
the error of the data points. We have estimated the cr from 



Fig. 1 of Alard (2001 1 as the dispersion of the star counts 



with respect to an exponential vertical profile. Afterwards 
we have used an iterative process to minimize using 
the Levenberg-Marquardt algorithm ( Press et al.||1992" l. 
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4.1. Fitting the disc and the boxy bulge 

We have not attempted to fit all the free parameters at 
once. First, we have fitted a star counts model with the 
disc and the boxy bulge to a reduced data set consisting 
of one point every 10 pixels (one data point every 0.5°). 
In the disc region, we have taken one point every 5 pix- 
els (0.25°) across the disc. Finally, we have neglected the 
data points in the inner region (—3 < l{deg) < 3, —2 < 
b{deg) < 2) since they are clearly dominated by the nu- 
clear bar. This reduced data set contains 1494 points and 
it is shown in Fig. |4] The errors of these data points used 
to calculate has been estimated to be 15 % as explained 
above. 

We have fitted the data for fixed values of the bulge 
angle (a^ = 60,45,35,20,10,0 deg) and the disc ra- 
dial scale {hr — 1.5,2,2.5,3.5 kpc). The results obtained 
for the other free parameters are shown in Table [l] For 
all a^'s, the best fits are always obtained with hr—2.5 
kpc. This is in agreement with the radial scale length 
of the disc derived from previous studies in the infrared, 
which have found values in the range from 1.9 (2MASS 
Lopez-Corredoira et al.|2005 ) to 2.6 (COBE Freudenreich 
1998||. Regarding the angle of the bulge, the best fit is 



obtained for 10°, although is only a 2.8 % higher for 
= 20° (Figure [i] shows the best fit with = 20°). 
This is also in agreement with other studies of the struc- 



ture of the inner Galaxy. Dwek et al. ( 1995 1 derived an- 
gles in the range 50° — 74° depending on the function 
assumed to represent the bulge (they found 74° with the 



same boxy Gaussian that we use here) , while Freudenreich 
( |1998 1 derived angles in the range of 75° — 81° also using 
COBE/DIRBE data. A high angle of 78° has also been 
measured from the 2 Microns Galaxy Survey (TMGS) by 
Lopez-Corredoira et al. (20001. On the other hand, the 
work of Lopez-Corredoira et al. (20051 using 2MASS data 
found a" = 20° -35°. 

The shape of the bulge as given by the best fit with 
= 10° is xq/xq : uo/xq : zq/xq = 1 : 0.5 : 0.3 while the 
best fit with — 20° has axis ratios of 1 : 0.55 : 0.4 
which are in perfect agreement with those derived by 



Lopez-Corredoira et al. ( 2005 ) also with 2MASS data or 



by Lopez-Corredoira et al. (20001 with TMGS data. In 



contrast, other studies have favored somewhat thinner 
bulges with ratios of 1 : 0.3 - 0.4 : 0.3 (|Dwek et al.||1995 



Freudenreich||1998| [Bissantz et al.||2003D 



4.2. Fitting disc, bulge and nuclear bar 

In the second data set we have also included the star 
counts in the innermost region of the Galaxy (—3 < / < 
3, —2 < 6 < 2). We have taken one point every 3 pixels 
(0.15°) across the disc. The data set is shown in Fig.js] It 
contains 1734 points. From Fig. 1 of Alard (2001 1 is clear 



that the dispersion of the star counts map in the nuclear 
region is smaller than in the bulge. We have estimated a 
a of 5% for these data points. 




-5 - 



10 







l(deg) 



-5 



Fig. 4. Thin solid contours represent the star counts map. 
Thick dashed lines are the best fit obtained with a star 
counts model with = 20° and hr = 2.5 kpc. The points 
show the data set used to fit disc and bulge. 



We have fitted the data using our star counts model 
with three components: disc, bulge and nuclear bar. We 
have taken as starting point the best fits for a^=10°, 20°, 
35°, 45° and 60° derived in the previous section (all the 
parameters of the disc and the bulge are fixed except the 
factors po). The best fit results for the parameters of the 
nuclear bar are summarized in Table |2] The full set of 
results for the nuclear bar fits are listed in Tables |3] - [H] 

The best fits are found with a thick triaxial bar with 
typical axis ratios of 1 : (0.7 — 0.8) : (0.5 — 0.6) and in 
the interval 60 — 120 for all the 's. This implies that the 
parameters of the nuclear bar can indeed be constrained 
independently of the bulge and that ^ . Globally, 
the best fit is obtained with = 10° and a nuclear bar 
with = 90° — 105°. some examples of good fits are 
shown in Fig. [5] 

Finally, we have also performed fits with a lopsided 
nuclear bar {xi ^ 0). The results are shown in Tables |3]-|6] 
and they are summarized in Table [2] The best fits are still 
obtained with a nuclear bar that is almost perpendicular 
to the Sun-GC line. The symmetry center of the bar would 
be shifted by 22-25 pc towards the third Galactic quadrant 
if a'' = 60° - 90° while it would be shifted by the same 
quantity but towards the fourth quadrant if = 90° — 
120°. 



5. Simulations of the gas dynamics: methods 



We have used the code of Combes & Gerin ( 1985 1 to sim- 



ulate the motion of gas clouds in the potential obtained 
from the mass distribution inferred with the star counts 
model. This code has already been used successfully to un- 
derstand the gas distribution and dynamics of a number 
of gal axies as MlOO (|Garcia-Burillo et al.||1998[) or NGC 
4736 dCerin et al.||l991| [Mulder fc Combes||1996 1. 
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Table 1. Results of the fits to the star counts data of Fig. [i] using an exponential disc and a triaxial boxy-Gaussian 
bulge 







Bulg 


e 






Disk 








yo 


zo 


Po 








Po 




kpc 


kpc 


kpc 




deg 


kpc 


kpc 






1.068 


0.498 


0.346 


0.389E-02 


0.0 


2.500 


0.168 


0.249E-03 


0.412 


1.045 


0.508 


0.348 


0.452E-02 


10.0 


3.000 


0.142 


0.365E-03 


0.290 


1.117 


0.513 


0.362 


0.383E-02 


10.0 


2.500 


0.153 


0.254E-03 


0.259 


1.243 


0.540 


0.383 


0.303E-02 


10.0 


2.000 


0.161 


0.145E-03 


0.264 


0.817 


0.491 


0.349 


0.615E-02 


20.0 


3.500 


0.134 


0.469E-03 


0.336 


0.840 


0.486 


0.359 


0.557E-02 


20.0 


3.000 


0.144 


0.366E-03 


0.288 


0.884 


0.489 


0.375 


0.477E-02 


20.0 


2.500 


0.156 


0.255E-03 


0.266 


0.975 


0.513 


0.400 


0.376E-02 


20.0 


2.000 


0.166 


0.145E-03 


0.284 


1.118 


0.574 


0.432 


0.283E-02 


20.0 


1.500 


0.169 


0.551E-04 


0.325 


0.703 


0.444 


0.355 


0.768E-02 


35.0 


3.500 


0.135 


0.470E-03 


0.335 


0.716 


0.435 


0.366 


0.708E-02 


35.0 


3.000 


0.145 


0.366E-03 


0.291 


0.744 


0.429 


0.382 


0.621E-02 


35.0 


2.500 


0.158 


0.255E-03 


0.274 


0.812 


0.442 


0.408 


0.499E-02 


35.0 


2.000 


0.169 


0.145E-03 


0.301 


0.937 


0.496 


0.445 


0.371E-02 


35.0 


1.500 


0.172 


0.551E-04 


0.350 


0.670 


0.381 


0.369 


0.851E-02 


45.0 


3.000 


0.146 


0.367E-03 


0.293 


0.693 


0.367 


0.385 


0.767E-02 


45.0 


2.500 


0.159 


0.255E-03 


0.278 


0.752 


0.365 


0.412 


0.637E-02 


45.0 


2.000 


0.171 


0.145E-03 


0.309 


0.621 


0.200 


0.373 


0.171E-01 


60.0 


3.000 


0.147 


0.367E-03 


0.297 


0.637 


0.132 


0.389 


0.226E-01 


60.0 


2.500 


0.161 


0.254E-03 


0.286 


0.681 


0.179 


0.416 


0.140E-01 


60.0 


2.000 


0.172 


0.144E-03 


0.327 



Table 2. Summary of the parameters of the nuclear bar as derived from the best fits to the 2MASS star counts data. 
The first row are fits with xi fixed to and the second row corresponds to fits with free xi to allow for lopsidedness. 
The full set of results as a function of the orientation of the bulge are listed in Tables [s] to Ib] 



Xq 




4 




a' 


Xi 




kpc 


kpc 


kpc 




deg 


kpc 




0.141-0.150 


0.130-0.138 


0.095-0.097 


0.141-0.150 


60-120 


0.0 


0.357-0.376 


0.158-0.166 


0.115-0.133 


0.088-0.094 


0.151-0.178 


60-120 


0.022-0.025 


0.331-0.346 



5.1. Stellar potential, rotation curve and masses 

The potential is obtained from the mass distribution de- 
rived in the previous section, using standard FFT tech- 
niques in a Cartesian grid of 512 x 512 x 16 cells. The 
cells are cubic with a linear size of 58.6 pc, therefore the 
physical size of the grid is 30 x 30 x 0.94 kpc'^. This is 
good compromise to model the large scale dynamics at 
least to a radius of 10 kpc and to have a good enough 
spatial resolution to study the Galactic center. 



As described by Combes & Sanders (1981), the inter- 



action between periodically reproduced images is avoided 
following the scheme of Hohl & Hockney (19691, i.e. the 



potential computations are done in a grid 8 times big- 
ger, of 1024 X 1024 X 32 cells. The disk mass is truncated 

at 15kpc in radius, and the mass distribution is assumed = f {p^ + + p^) + p 
spherical outside, so that its influence in the center is 
cancelled. This is equivalent to assume that a spherical 
dark halo is dominating the mass outside. This assump- 
tion should have a negligible influence on the dynamics of 
the inner Galaxy, which is the main interest of this paper. 



The stellar potential is decomposed in its axisymmetric 
and non-axisymmetric parts. The axisymmetric potential 
is calculated in radial bins using the averaged mass in 
the bin. The non-axisymmetric part of the potential is 
calculated with a mass distribution corresponding to the 
difference between the total-mass and the axisymmetric 
distribution. 

The stellar number density derived in the previous sec- 
tion should be multiplied by a factor / to obtain the total 
mass density. This factor accounts for the mass of the stars 
detected by 2MASS but also for the mass of those stars 
that are not counted in the 2MASS maps. In addition, 
we have to introduce a dark halo to explain the observed 
velocities at large radii. The total density used is then: 



H 



(18) 



where the halo density is defined by a Plummer law: 



[1 -f- (r/r,)2]5/2 



(19) 
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Fig. 5. Upper panels: Thin sohd contours represent the star counts map. Thick dashed hues are the best fit obtained 
with a star counts model with = 10°, hr = 2.5 kpc and = 75 deg. The panel on the right is a zoom of the 
nuclear region. The points indicate the data used to do the fit. Lower panels: Same that the upper panels for the best 
fit with = 45°, and a lopsided bar with a'' = 60°. 



To determine the scaling factors / and Pq we have 
compared the terminal velocities and the rotation curve 
of the models with those derived from spectroscopic ob- 



servations (Clemens 1985 Fich et al. 1989 Burton & 



Liszt 19931. To obtain terminal velocities of the right 
magnitude (see Figs, [s] and |9| and a flat rotation curve 
at large radii (Fig. [6f7 one should apply a scaling fac- 
tor / = 3.9 10* and a halo with = 1.1210^ and 
Tc — 15.0 kpc. With these factors and the parameters 
given in Tables [3] - [g] the function gives the total 
mass density in units of 10^ MQ/kpc"^. Integrating in the 
total grid size, the typical masses of the different com- 
ponents are: = 710^° M©, = 2.9 10^° M©, 
= 1.9101" M© and = 2.0 10^ M©, which are 
in good agreement with previous determinations (see for 
instance [Weiner fc Sellwood||1999[ |Fux||1999[ ). 



Figure |6]shows an example of rotation curve. The mod- 
eled curve compares well with the rotation curve derived 
from observations of the interstellar gas for radius larger 
than 2 kpc. The agreement seems worse for small radii. 
However the comparison in this region is not straightfor- 
ward. On the one hand, the measured curve has been de- 
rived using the tangent point method assuming that the 
clouds move in circular orbits, which is not true in the 
innermost region. In addition, this curve depends on the 
position of the observer in the Galaxy. On the other hand, 



the modeled rotation curve has been derived with the ax- 
isymmetrical potential and with an azimuthal average of 
the non-axisymmetrical potential. This rotation curve is 
a kind of average curve independent of the position of the 
observer. 

Figure [7] gives the rotation velocity as a function of 
galactocentric radius. It also shows the + n/m curves 
(where k is the epicyclic frequency). For a given bar pat- 
tern speed Hp, these curves give the location of the corota- 
tion {il = Hp), the Outer Lindblad resonance (il -I- k/2 = 
Hp), the Inner Lindblad resonances (H — k/2 = Hp) and 
the Ultra Harmonic Resonance {H — k/4 = Hp). 



5.2. Gas cloud dynamics 

The gas clouds are modeled as particles that are initially 
launched in the axisymmetric potential with a velocity 
corresponding to the circular velocity of the potential at 
the radius of the particle. The initial radial distribution 
is an exponential disk with radial scale of 2 kpc. The 
non-axisymmetric part of the potential is then introduced 
gradually multiplying by the factor 1 — exp[— f/r]^ with a 
time delay of r =100 Myr. The angular speed of the non- 
axisymmetrical pattern. Hp, is constant along the simula- 
tion. 
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Fig. 6. Rotation curve obtained with the halo and the 
scaling discussed in Sect. 5.1 for a bulge and nuclear bar 
orientation of = 20° and = 75°, respectively. 
The different lines show the contribution of the differ- 
ent mass components to the total rotation curve. The red 
points indicate the rotation curve measured from obser- 
vational data. The mass of the halo , disc, and bulge 
are M" = 710^° Mq, M° = 2.9 lOi"^ Mq, = 
1.9 10^" M(n, respectively. The mass of the nuclear bar 



is M" 



2.0 10^ M. 
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Fig. 7. Angular frequencies as a function of radius in the 
epicyclic approximation for the potential discussed in Sect. 
5.1 



& Gerin (1985) and Casoli & Combes (19821 



6. Simulations results and discussion: large scale 
structure 

Figures [8] and [9] show face-on views and Zw-diagram for a 
grid of models with from 10° to 45° and pattern speeds 
from 20 to 50 km s^^ kpc^^. For an easy comparison with 
the observations, we have overlaid in the Zv-diagrams the 
TVC points, the boxes indicating the tanget points and 
the solid and dashed lines shown in Fig. [l] We have also 
plotted the four logarithmic spirals of |Nakanishi fc Sofue| 
(2006) in the face-on views. The models shown in Figs. 
|8| and [9] exhibit a number of common characteristics. The 
galaxies have two major spiral arms but they also exhibit a 
second pair of less prominent spiral arms. In the inner ~ 4 
kpc, there are arms that contour the bar forming a kind 
of oval (hereafter lateral arms) . The density of gas inside 
this oval decreases down to a radius of ^ 1 kpc, where 
there is a ring. This ring is connected to the lateral arms 
by the off-axis arms (or dustlanes). In the inner hundreds 
of parsecs there is a nuclear condensation of gas. 

The Milky Way galaxy is a very complex system and 
up to now there is not a single simulation of the gas dy- 
namics that can explain all the observed features, nor 
quantitatively nor qualitatively. [Englmaier fc Gerhard] 
(1999) have compared their models to the terminal ve- 
locity curve, they have also compared the position the 
arms tangent points in the CO /w-diagram, at |^| — 30° 



and 50°. The two Zw-diagrams presented by Fux ( 1999 1 has 
been selected by the global resemblance to the ^w-diagrams 
obtained from CO and HI data. They reproduce qualita- 
tively the 3-kpc arm and the Connecting Arm. |Weiner| 



& Sellwood ( 1999 1 constrained their models by compar- 



ing with the HI terminal velocity curve. Finally, |Bissantz" 



et al. ( 2003 1 have used mainly the terminal velocity curve 
and some voids in the observed ?w-diagram to compare 
with their simulations. To select the best models in the 
present work we have compared with: i) the terminal ve- 
locity curve, ii) the inclination of the GMR in the Iv- 
diagram, in) the existence or not of the 3kpc arm and 
when it exists, its inclination and absolute velocity iv) the 
arms tangent points v) the shape and size of the HI Ikpc 
ring. 



We do not use a mass spectrum for the clouds and 
we do not simulate star formation. The simulations con- 
tain typically 10^ particles. The effects of the gas self- 
gravitation are neglected. The clouds move as test parti- 
cles in the stellar potential except when they collide in- 
elastically. The forces due to the stellar potential are eval- 
uated every 1 Myr and the particles are moved according 
to them. The particles are followed in three dimensions 
but the collisions are computed in a two-dimensional grid 
with 480 X 480 cells. Clouds in the same grid cell can col- 
lide. The relative velocity of two collisions partners loses 
60% of its absolute value in the coUision. A detailed dis- 



6.1. Pattern speed and orientation of the bulge 

Figure |8] shows face-on views and Iv-diagrams for a grid of 
models with — 20° and pattern speeds from 20 to 50 
km kpc~^. The model with Up — 20 km s~^ kpc^^ 
reproduces some spiral tangent points and has an inner 
lateral arm with the inclination of the 3-kpc arm, how- 
ever it does not reproduce the terminal velocity curve nor 
the velocity of the 3-kpc arm, for instance. On the other 
extremity, the model with fip = 50 km s~^ kpc~^ does 
not reproduce satisfactorily neither the terminal velocity 
curve nor the GMR nor the 3-kpc arm. 
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In contrast, the models with flp — 30 and 40 
km s^^ kpc~^ give a very good overall representation of 
the Galaxy. Both reproduce satisfactorily the spiral tan- 
gent points at \l\ = 30,50° within ±5° and ±10 km s~^. 



approximately. In addition, the model with ftp 



40 



km kpc~^ reproduces very well the Carina arm and 
its tangent point at « —80°. However, the locus of the 
inner spiral arms and the 3-kpc arm in the modeled Iv- 
diagram are too steep in comparison with the observed 
Zf-diagram. In contrast, the GMR and the the 3-kpc arm 
is very well reproduced from |/| = 90° to 15° by the model 
with flp = 30 km kpc~^. In particular the 3-kpc 
arm is reproduced with both the good inclination and the 
velocity within ^10 km s~^. In addition, the terminal 
velocity curve for / > 5° is very well reproduced. Given 
the difficulty of the task, it will not be realistic to give a 
very narrow interval for the bar pattern speed, from our 
best models we conclude that the most probable value for 
Op should be in the range 30 — 40 km s~^ kpc^^. However, 
since this work is devoted to the inner Galaxy structure 
and kinematics in the following we will mainly discuss 
models with Op — 30 km s~^ kpc~^ since they explain 
the terminal velocity curve, the 3-kpc arm and the GMR 
better than the model with Op = 40 km s~^ kpc~^. 

Figure |9] shows face-on views and Zu-diagrams for a 
grid of models with constant O = 30 km s~^ kpc"-'^ and 
a^'s of 45° , 35° , 20° and 10° {a^ is constant to 90° 
for all the models). For = 10°, the velocity along the 
line of sight is too low to explain the high non-circular 
velocities of the 3 kpc arm (—53 km s~^ at I — 0°) and 
the observed arm is not reproduced in the Iv diagram. In 
contrast, the 3-kpc arm and the inclination of the GMR 
is well reproduced for from 20° to 45°. Nevertheless, 
the model with = 45° does not reproduce correctly 
the positive velocity terminal curve. In addition, there is 
no ring at 1 kpc in the face on view. 

Once again, we do not pretend to give a very narrow in- 
terval of angles for the inclination angle of the bar since the 
analysis, as in other published works, is still rather qual- 
itative. We reckon that the inclination of the bar should 
be in the range from 2(F to 35° , with some preference for 
20° since the quantitative agreement of the 3-kpc arm and 
the terminal velocity curve is better. Therefore, the model 
with Op = 30 km s^^ kpc""'^ and — 20° will be con- 
sidered as our standard model in the rest of the paper. 



and stars simulations, Fux found that the pattern speed 
evolves from 50 to 30 km s^^ kpc~^. Taking into account 
our potential, a pattern speed of 30 or 40 km s~^ kpc~^ 
gives a corotation radius of 7 and 5 kpc, respectively (see 
Fig.[7|. Therefore, the corotation radius is well beyond the 
bar extremity, which is the characteristic of a slow bar (see 
for instance |Sellwood fc Sparke||1988 1 . 

To explain the fact that the dustlanes in the inner 
spiral arms of some galaxies are located in the concave 
section of the arms, and assuming that the bars are fast 



(i.e., they extend to the corotation), Sellwood & Sparke 



( 1988 ) proposed the existence of different speeds for the 



bar and the spiral patterns. Indeed, for the Milky Way 
some works have derived low rotations speeds of 13-20 
km s~^ kpc~^ for the spiral pattern (Lin et al. 1969 



Morgan 1990 Amaral & Lepine 19971. Gas flow models 



with independent speeds for the bar and the spiral pat- 



terns have been presented by Bissantz et al. (20031, who 



discussed that different speeds for the spiral and the bar 
patterns could explain some regions devoid of gas in the 
Zu-diagram. However, the overall agreement of their Iv- 
diagrams with the observations is not better than other 



models like that of Fux (19991 or our models, which are 



able to reproduce quantitatively features as the 3-kpc arm. 

It is interesting to remark that our simulations favor 
a relatively slow bar and not a fast bar as assumed by 



Sellwood & Sparke ( 1988 1 . Furthermore, the pattern speed 
of the relatively slow bars models (30-40 km s^^ kpc~^ 



Weiner fc SeUwood|[T999| |Fux||1999[ this work) IS compa- 
rable to several determinations of the spiral pattern speed 



that give values of 30-35 km s kpc (Mishurov & 



Zenina,1999,, Fernandez et al.|2001 [Lepine et al.|2001[ ). On 



the other hand, Ibata & Gilmorel ( 1995 1 have measured a 



rotation speed of 25 ± 4 km kpc~^ for the bulge in 
the 0.7 < R < 3.5 kpc region. Therefore, at present there 
are still many uncertainties on the speeds of the differ- 
ent patterns and it is not clear whether the spirals and 
the bar are decoupled or they turn at the same speeds. 
Simulations with an independent pattern speed for the 
spirals would be needed to study this question in detail 
but this is beyond the scope of the current paper, which 
is mainly devoted to the inner Galaxy. 



6.2. The GMR 



6.1.1. Comparison with previous works 

The inclination of the bar derived from our models is 



in good agreement with previous determinations (Binney 



[eTal] [19971 [FiEl [19991 |L6pez-Corredoira et al.|[2005 



see 



also Sect. 1). Regarding the bar pattern speed, our simu- 



lations give results that are lower than those of Englmaier 



fc Gerhard (19991 or Bissantz et al. (20031, which sug 



gest Op = 50 — 60 km s ^ kpc ^, but in good agreement 
with the results o f [Weiner fc Sellwood] ( [1999| , Op = 42 
km s^^ kpc^^, or Fux (19991. In his self-consistent gas 



To get further insight in the structure and the dynamics 
of the Galaxy, we present below a more detailed analysis 
of the model with Op = 30 km s~^ kpc~^ and a^=20°. 
We have selected different structures in the face-on view 
and we show the locus of the different structures in the 



Zii-diagram (Fig. 10 1 



As already mentioned in Sect. 2.2, it is uncertain 
whether the GMR is an actual ring or it is composed by 
imbricated spiral arms. In our models the radius of the 
GMR (4-5 kpc) corresponds exactly with the position of 
the Ultra-Harmonic resonance for O ~ 30 km s~^ kpc~^ 
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(Fig. I?]). However, we do not support the idea that the 
GMR is a resonant ring. 

As shown in the upper panels of Fig. 



10 in our mod- 



els the GMR is composed by the inner parts of the spiral 
arms, which at those radius show almost circular motion. 
This is in agreement with the face-on views of the CO and 



HI data of ,Nakanishi fc Sofue ( 2006, ), who have proposed 
that the GMR is formed by the Scutum-Crux arm (near 
side of the GMR) and the inner part of the Sgr-Carina (far 
side of the GMR) . .Englmaier fc Gerhard] ([l999 ) discussed 
that if formed by imbricated arms, the GMR should be 
outside the corotation assuming that there cannot exist 
significant circular orbits inside the corotation to produce 
a feature resembling a circular ring. However, this state- 
ment only holds for fast bars, for which all the region 
inside the corotation is dominated by the bar. In contrast, 
in our models with a relatively slow bar, the inner part of 
the spiral arms is inside the corotation radius and exhibit 
almost circular motions. 

6.3. The 3-kpc arm and its far side counterpart 

Lower panels of Fig. [TO] show the locus of the lateral arms 
that contour the bar. The two dot-dashed lines in the Iv- 
diagram are the fit to the 3-kpc arm and to its far-side 
counterpart as given by Dame fc Thaddeus (20081. As 



already mentioned, the locus of the lateral arms in the 
/i;-diagram reproduce quite well the structure of the 3- 



kpc arm. This is in agreement with Fux (19991, who has 



already proposed that the 3-kpc arm could be a lateral 



arm. On the other hand, Fux ( 1999 1 proposed that the 
"135 km s~^ arm" is the far side counterpart of the 3-kpc 
arm. However, our models predict a far side counterpart 
of the 3-kpc arm whose locus in the Zii-diagram is almost 
symmetrical to the 3-kpc arm. This is in perfect agreement 



with the recently found far 3-kpc arm (Dame fc Thaddeus 



20081. The full implications of the new finding by Dame 



fc Thaddeus (20081 will be presented elsewhere. 



The fact that the 3-kpc arms are lateral arms im- 
ply that they cannot be a channel to transport gas from 
the corotation to the CMZ as suggested by [Habing et al.| 
( 2006 1 , at least not directly. 



6.4. The off-axis arms or dustlanes 

In contrast to the 3-kpc arms, the dustlanes do contribute 
to the transport of gas to the inner regions. In our mod- 
els the locus of the dustlanes in the /w-diagram are verti- 
cal structures with velocities from ^ km s~^ to more 
than 100 km s~^. The maximum velocity increases as I 
decreases. These vertical structures in the ?w-diagram re- 



semble the Clumps found by Bania (19771 and discussed 
recently by ,Liszt] ( 2006[ 20081. Nevertheless, our models 
do not reproduce the Galactic longitude of the observed 
Clumps (5.5° and 3.1°), we confirm that Banias's Clumps 
are probably shocked gas in the dustlanes, which was the 



"Clump 2" . The fact that the Clumps are composed of 
gas suffering the dustlanes shocks is also suggested by ob- 
servational evidences both kinematical (linewidths, veloc- 
ity gradients, Stark fc Bania||1986 Liszt||2006 ) and chem- 



ical ( Rodriguez- Fernandez et al. 



20061. Taking into ac- 



count the discrete appearance of the Bania's clumps at 
two Galactic longitudes, it is likely that the dust lanes of 
the Milky Way are patchy as it is commonly observed in 
external galaxies. 



Fux (19991 and Marshall et al. (2008) have proposed 



a slightly different interpretation. They propose that the 
Connecting Arm is the locus of the near side dustlane 
in the Zu-diagram. However, the linewidths of the gas in 
the Connecting Arm do not show the velocity disper- 
sion expected in the dustlanes shocks (see for instance 
Reynaud fc Downes 1998) . Indeed probably both the 



Bania's clumps and the Connecting Arm are related to 
the dustlane. We reckon that Bania's Clumps are clouds 
that are actually suffering the strong shocks expected in 
the dustlanes while the Connecting Arm is most likely 
the post-shocked and accelerated gas (see also the discus- 
sion in [Liszt 2006) . Indeed, in our simulations the max- 



imum velocity of the gas in the dustlanes increases as I 
decreases with the same steep of the Connecting Arm. 
Observationally, this interpretation is supported by the 
fact that the Clump at I — 5.5° seems to be connected to 



to Connecting Arm at a negative latitude (see Bitran et al. 
1997). The situation was less clear for the "Clump 2" at 
I = 3.1° since it is found at positive latitudes. However, 
this apparent problem has been solved with the recent 



discovery by Liszt ( 2008 1 of an emission feature with the 



same inclination in the Zu-diagram that the Connecting 
Arm but at positive latitude. In addition it seems con- 
nected to the Clump 2. In any case, the Clump 2 exhibits 
a rich shock chemistry as expected for the gas in the dust- 
lanes shocks ( Rodriguez-Fernandez et ar]|2006 l 



original interpretation of Stark fc Bania (1986) for the 



7. Simulations results and discussion: nuclear bar, 
HI ring and CMZ 

Figure [Tl] shows the face-on view of the inner 2 kpc of the 
modeled galaxies and the Zw-diagrams for different orien- 
tations of the nuclear bar. All the other parameters are 
the same that those discussed in the previous section and 
the large scale face on views and Zw-diagrams are the same 
independently of the orientation of the nuclear bar. 

The dustlanes end in a ring and inside this ring there 
is a small bar-like structure that is connected to the ring 
by two small spiral arms. The HI ring has inner and outer 
radius of 300 pc and 800 pc, respectively, with are in excel- 
lent agreement with the HI observations of |Liszt fc Burtori] 
(1980) (see their figure 3). The ring is almost circular. In 
some simulations we have also found elliptical rings with 
the major axis approximately perpendicular to the large 
bar. Therefore, the ring is not supported by xl orbits elon- 
gated in the direction of the large bar (see a review of the 
different proposed models in Ferriere et al.||2007 ). 
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7.1. The orientation of the nuclear bar 

In the central hundreds of parsecs there is an elongated 
structure that corresponds to the gas response to the nu- 
clear bar and that resembles the observed CMZ. 

We have analyzed the shape of the Zw-diagram as a 
function of the nuclear bar orientation. Figure [TT] shows 
the Zi!-diagram obtained for nuclear bar orientations from 
0°, to 15°, 30° and 45°. Our modeled Zu-diagrams do 
not reproduce the components K and J of Rodriguez- 
Fernandez et al. (2006). In contrast, the Zw-diagram of 
the CMZ is a kind of parallelogram that resembles the ob- 
served one (Fig. |2|. The inclination of the observed par- 
allelogram is very well reproduced for an angle of 15°. 
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B'™ey| ( |1994| tried to reproduce the parallelogram dis- pjg_ ^2. Same that Fig. g but with a nuclear bar mass of 



cussed by Binney et al. (1991) using n-body simulations 



but with little success. To our knowledge, this is the first 
time that numerical simulations of the gas dynamics of 
the Milky Way reproduce naturally the parallelogram of 
the ?u-diagram of the CMZ. This is probably due to the 
accurate potential that we have determined with the deep 
2MASS data. 

It is worth noting that the orientation angle of the 
small bar inferred from our simulations is in excellent 
agreement with our models of the 2MASS star counts and 
with the face on view of the CMZ inferred from CO and 
OH data by ISawada et al.l (|2004|). We conclude that the 



observed CMZ is most likely the gas response to the nu- 
clear bar and that the orientation of the nuclear bar is 
- 75°. 



7.2. The mass of the nuclear bar 

In the simulations discussed above, the velocity dispersion 
of the CMZ is well reproduced for a mass of the nuclear 
bar of 2.010^ M© (Sect. 5.1). However, the terminal ve- 
locities of the HI ring do not match the observed ones. 
To get further insight on the nuclear bar mass we have 
increased its value in order to reproduce the terminal ve- 
locities of the HI ring. In these models, the shape of the 
different mass components is still given by the fits to the 
2MASS star counts map (Tables [3]-|6|. However, we have 
considered different scaling factors for the nuclear bar and 
for the bulge and disc. Thus, instead of Eq. [I8] the total 
density is given by: 
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(20) 

Therefore, with respect 



where p is given by Eq. 
to previous models, the only parameter that we change is 
/''. Figure 12 shows the rotation curve for f'' = 1.1 10^, 
which gives a nuclear bar mass of 5.510® Mq . All the other 
parameters and masses are the same that those mentioned 
in Sect. 5.1. The simulations results for the inner 2 kpc 
of the Galaxy are shown in Fig. [13] The face on view is 
characterized by a large ring with inner and outer radius of 
200 and 1200 pc, respectively. The Zt;-diagram of this ring 
reproduces the terminal velocity of the HI observations. 



M'' = 5.5 10® Mr. 



which at positive longitude coincides with the Connecting 



Arm. Indeed, Liszt & Burton ( |1980) proposed that the 
Connecting Arm could be the edge of the HI ring. 

It is difficult to reproduce both the terminal velocities 
of CMZ and HI ring, however from the previous considera- 
tions we conclude that the mass of the nuclear bar should 
be in the range (2 - 5.5) 10® Mq. 

7.3. Is the nuclear bar lopsided? 

The nuclear bar in the star counts map of Alard| (2001) 
seems to be shifted towards negative longitudes (in con- 
trast to the gas distribution that seems to be shifted to- 
wards positive longitudes). Our star counts model shows 
that fits with a lopsided nuclear bar have a that is 
9% lower than the for models that do not allow a lop- 
sided bar. The symmetry center of the bar could be shifted 
by 22-25 pc with respect to the dynamical center of the 
Galaxy. 

We have studied the gas response to a lopsided nuclear 
bar in order to investigate whether the observed asymme- 
try of the CMZ could be explained by a lopsided nuclear 
bar as inferred from the 2MASS data. As discussed by 



Morris & Serabyn ( 1996 ) the signature of a m = 1 mode 
in the Galactic center would be a shift of the gas distri- 
bution in I but also a shift of the velocity centroid of the 
ZiJ-diagram. 



Figure 14 shows a simulation of the gas flow in the po- 
tential computed from the flt to the 2MASS star counts 
map with a* = 75° (for which xi = 0.023 kpc). The Iv- 
diagram of the gas exhibit a parallelogram whose sides 
show the good inclination. The simulated parallelogram 
is only slightly asymmetric in I and v. To better under- 
stand the effect of a lopsided nuclear bar, we have also 
computed a simulation with the same nuclear bar but ar- 
tificially shifted by 90 pc along its major axis towards the 
third Galactic quadrant {xi = 0.090 kpc). The results are 
also shown in Fig. [14] The Zv-diagram is still only slightly 
asymmetric in v but now it is clearly asymmetric in I . The 
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gas response to the lopsided stellar bar follows the star dis- 
tribution and it is lopsided towards negative longitudes as 
well. To obtain a ^w-diagram with a parallelogram sim- 
ilar to that of the CO data shown in Fig. [2] the nuclear 
bar should be lopsided towards the first quadrant. This is 
shown in Fig. 14 where xi = —0.090 kpc. However, such 
a bar is not compatible with the 2MASS star counts. 

We conclude that, both from the fits to the 2MASS 
data and from the numerical simulations of the gas dy- 
namics, there are no clear evidences of an intrinsic lopsid- 
edness in the stellar potential. The observed asymmetry 
of the gas distribution in the CMZ cannot be explained 
as the gas response to a lopsided nuclear bar in the way 
suggested by Alard's map. 



8. Discussion 

8.1. The origin of the asymmetry of the CMZ 

Although the origin of the asymmetry of the CMZ is a 
long-standing problem, it has not been discussed much in 
the literature. Two exceptions are the papers by |Combes| 
(19961 and Morris & Serabyn (19961. The first paper pro- 



posed that the observed asymmetry could well be due to 
a, m — I mode of the potential since lopsidedness is com- 
monly observed in the central regions of external galaxies. 
The second paper discusses how the effect of the pres- 
ence of a m = 1 mode in the Galactic center would be an 
asymmetry of the CMZ both in longitude and in velocity. 
However, we have shown in the previous section that the 
possible lopsidedness of the stellar distribution cannot ex- 
plain the observed asymmetry of the CMZ. Furthermore, 
there is not a clear evidence of an intrinsic asymmetry of 
the stellar potential. Regarding the gas, the distribution 
is clearly asymmetric in I but the parallelogram of Fig. 
[2] is rather symmetric in v. In addition, the structures 



found by Sofue ( 1995 ) that seem to be the signature of a 



quasi-circular ring inside the CMZ and that could be the 
nuclear bar equivalent to the HI ring, do not show neither 
an asymmetry in v. 

Indeed our simulations propose an alternative explana- 
tion for the asymmetric distribution of the CMZ in which 
the stellar potential in the GC is strictly symmetric. The 
simulations show that the HI ring can be connected to 
the CMZ by a pair of small spiral arms. For instance, the 
simulation with a* = 75 (second row of Fig. 11) shows 
that the clouds in the arm seen at negative longitudes will 
be seen in the /u-diagram outside of the observed paral- 
lelogram while the clouds in the arm seen at positive lon- 
gitudes "fill" the observed parallelogram at the position 
of the 1=1.3° -complex. Observations of the molecular gas 
in external galaxies show that, in contrast with most nu- 
merical simulations, the gas distribution in real galaxies is 
rarely symmetric. This could also be the case in the Milky 
Way. The apparent asymmetry of the CMZ can be due to 
material falling into the CMZ from the HI ring through 
only one of the inner spiral arms. In this context, strong 



shocks are expected in the interaction region, that would 
be the 1=1.3^ -complex. 

Observationally, the 1=1.3!^ -complex looks very differ- 
ent to the Sgr's cloud complexes. For instance, in the Ih 
map of Fig. [2] it is clear that the gas distribution from Sgr 
C to Sgr B is rather symmetric and thin. The asymme- 
try in I is indeed exclusively due to the 1=1.3° -complex. 
Furthermore, the 1=1.3° -complex shows a very special lat- 
itude extension, much higher that the rest of the CMZ. 
The velocity structure of this cloud complex is also quite 
singular, showing a large velocity dispersion (the 1=1.3° - 
complex is indeed the most prominent feature in the Iv- 
diagram). The chemistry gives us also interesting hints 
on the nature of this cloud complex. On the one hand. 



Hiittemeister et al. (19981 have found that the highest 



abundances of the shock tracer molecule SiO, are found at 
one of the extremities of the CMZ: in the 1=1.^?^ -complex. 
They have already proposed that this cloud complex is the 
scenario of strong shocks due to gas falling into the CMZ. 



On the other hand, Rodriguez- Fernandez et al. ( 2006 1 



have presented a more complete study of the SiO emis- 
sion in all the kinematical structures of the GC. They 
have detected SiO in at least one cloud of every kinemat- 
ical structure. However, the only feature where they de- 
tected high abundances of SiO in every observed cloud is 



the 1=1.3° -complex. In agreement with Hiittemeister et al. 
(1998), Rodriguez-Fernandez et al. (2006) have measured 



the highest SiO abundances in this region. Therefore, the 
lb distribution, the velocity structure, and the chemistry, 
imply that the cloud complex giving rise to the observed 
asymmetry of the CMZ (the 1=1.3° -complex) shows signif- 
icant differences to other CMZ complexes. In particular, it 
shows the signature of strong shocks. Therefore, we reckon 
that the observed asymmetry of the CMZ can well be the 
result of gas falling into the CMZ and not due to an intrin- 
sic asymmetry of the stellar potential. This explanation is 
also supported by our detailed simulations of the stellar 
structure and the gas dynamics. 

8.2. The nature of the nuclear bar 

Double bars are commonly observed in external galaxies 



1996 



(Shaw et al. 1993 Wozniak et al. 1995 Elmegreen et al 



Around 28 % o f barred galaxies have a secondary 
bar ( |Erwin fc Sparke||2002| |Laine et al.]|2002[ ). Nuclear or 
secondary bars are objects of prime interest to explain the 
gas inflow towards the center of a galaxy and the fueling 
of Active Galactic Nuclei ( Shlosman et al.||l989 l. 

Several works have studied the stellar orbits that sup- 



port nested bars (Maciejewski & Sparke 1997 2000) and 



their stability ( El-Zant & Shlosman 2003 1 . In addition 



double barred galaxies have been the subject of a number 
of numerical simulations of the star and gas dynamics that 
have shown a plethora of phenomena as: i) secondary bars 
forming in purely stellar discs ( Rautiainen et alT]|2002 1 or 
only in the presence of gas ( Combes||1994 ). ii) secondary 
bars forming before ( Rautiainen fc Salo|1999[ ) or after the 
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primary bars ( 


Combes|1994 Englmaier & Shlosman|2004 


Heller et al.||2007l, Hi) decoupled secondary bars rotating 


slower (Heller et al. 


20011 or faster (|FriedH & Martinet 


1993 Combes||1994| 


■lautiainen et al.||2002| |Englmaier & 


Shlosman|2004 


1 than the primary bar, iv) two misaligned 



19941 



The secondary bars rotating at the same speed than 
the primary bars have been explained by |Shaw et al. 
(19931 and Combes (19941 as follows: the gas clouds tend 



to shift from the closed xl orbits to the perpendicular 
x2 orbits of the primary bar. However, dissipative colli- 
sions between the gas clouds reduce the orthogonality of 
this phase shift and the gas settles in a leading phase- 
shifted bar. If the gas fraction is high, its gravitational 
influence is then sufficient to modify the stellar compo- 
nent and to form a secondary bar, that is still coupled to 
the primary bar. If secondary and primary bar remain cou- 
pled and they continue rotating with the same speed, the 
bars will align or form an angle of 90° since the primary 
bar exerts a gravitational torque on the secondary bar 



and vice- versa (Friedli & Martinet 19931. Alternatively. 



the secondary bar may increase its rotation speed. Thus, 
the coupled phase can be a precursor of a future decoupled 
phase. The dynamical process from coupled to decoupled 
gaseous nuclear bars have been studied by [Heller et al.| 
(20011, Englmaier & Shlosman (20041 and, more recently. 



by Heller et al. ( 2007 1, who have done the first simulations 
of the formation of nested bars from cosmological initial 
conditions. In these simulations a secondary bar forms in 
response to the gas inflow along the stellar primary bar. 
This gaseous bar is initially corotating with the primary 
bar. Subsequent mass inflow strengthens this bar giving 
rise to even more rapid gas inflow. Along this process the 
secondary bar shrinks size and the pattern speed increases 
inversely proportional to the bar size. The pattern speeds 
of the two bars are such that the secondary bar corotation 
coincides with the primary bar Inner Lindblad Resonance, 
suggesting that non-linear interactions between the two 
bars are at play ( Tagger et al.||1987 l. 



Our simulations are the very flrst attempt to model the 
Milky Way with two nested bars and have been done as- 
suming a common pattern speed for both bars. Therefore, 
we cannot study in detail the coupling and the evolution 
of the nuclear bar. Nevertheless, it is worth noting that 
these simulations explain many characteristics of Galaxy 
at scales from the disk to the nuclear region, some for the 
first time as the parallelogram of the CMZ. Therefore, the 
observational data are compatible with a scenario of cou- 
pled bars rotating with the same speed. Our results imply 
that the nuclear bar is leading the large primary bar by 
~ 55°, as expected in this context ( Combes||1994 l. Thus, 
the decoupling of the nuclear bar can still not be effective 
in the Milky Way. The dynamical decoupling of the two 
bars will be the subject of a forthcoming paper. 



9. Conclusions 

We have presented gas flow models in the mass distribu- 
tion derived from 2MASS star counts using a model with 
three components (disk, bulge and nuclear bar). Our dy- 
namical models are the first ones that include a central 
mass cusp (the nuclear bar) constrained by observations. 
For the first time, we have obtained good models of the 
Milky Way from the scales of the spiral arms to the Central 
Molecular Zone (CMZ) in the Galactic center (GC). The 
best models are found for a bulge orientation of 20 — 35° 
with respect to the Sun-GC line and a pattern speed of 
30-40 km s^^ kpc^^. Which places the corotation radius 
at 5-7 kpc. The simulations reproduce: 

— the spiral arms, giving in particular the good tangent 
point for the Carina arm. 

— the Galactic Molecular Ring (GMR), which is not an 
actual ring but the inner parts of the spiral arms and 
it is found approximately at the position of the ultra- 
harmonic resonance 

— the 3-kpc arm and its far-side counterpart, which are 
the lateral arms that surrounds the bar 

— the HI ring, with a very similar size to that derived 
from the observations 

— the Central Molecular Zone (CMZ) 

In addition, we postulate the existence of small spiral arms 
arising from the extremities of the CMZ and linking the 
HI ring. 

In our simulations, the CMZ is the gas response to 
the nuclear bar. These simulations reproduce, for the first 
time, the parallelogram shape of the Zw-diagram of the 
CMZ. Using this shape we have been able to constrain 
the nuclear bar orientation, which is of ~ 60° — 75° with 
respect to the Sun-GC line, in excellent agreement with 



the 2MASS fits and the results of [Sawada et al.| ( |2004 |. 
We have also studied the observed asymmetry of the CMZ 
and we have concluded that it cannot be due to the lopsid- 
edness of the nuclear bar suggested by the 2MASS maps. 
Gas dynamical simulations in a lopsided potential cannot 
reproduce the observations, and in any case, we do not 
see a strong direct evidence in 2MASS data of an intrinsic 
lopsidedness of the stellar potential. 

We propose an alternative scenario to explain the ob- 
served asymmetry of the CMZ. The asymmetry can be 
due to gas falling into the CMZ trough the inner spiral 
arm in the 1=1.3^ -complex region. We have also discussed 
that the 1=1 .3° -complex shows significant differences with 
the, otherwise symmetric, CMZ that extends from Sgr C 
to Sgr B. In particular, the 1=1. 3° -complex shows all the 
signatures of the shocks expected if it is an interaction 
region where new material is falling into the CMZ. 

In our models both bars rotate with the same pattern 
speed. The success of these models in explaining many 
characteristics of the Milky Way as the Central Molecular 
Zone (CMZ) imply that the observations are compatible 
with that assumption. Furthermore, in our models the nu- 
clear bar is leading the large bar by ~ 55° as expected for 
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coupled nested bars ( Combes|1994 l. Therefore, in spite of 



the high mass of the nuclear bar (10-25 % of the mass of 
the bulge) the dynamical decoupling my still not be effec- 
tive. However, this question requires further studies and 
gas flow models with two pattern speeds or self-consistent 
gas and stars dynamics simulations. 
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Fig. 8. Simulations results for — 20°. The different rows represent the face-on view and the longitude velocity 
diagram for different pattern speeds of 20, 30, 40 and 50 km/s/kpc (from upper to lower rows). The solid lines represent 
the spiral arms as defined by 'Nakanishi fc Sofue| p006| . The hne Sun-GC first crosses the Sgr-Carina arm and the 
Scutum-Crux arm. On the other direction, the line of sight to the Galactic anticenter crosses the Perseus arm. The 
fourth arm is the Norma-Outer arm. 
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Fig. 9. Simulations results for Hp = 30°. The different rows represent the face-on view and the longitude velocity 
diagram for different orientations of the bulge of 45, 35, 20 and 10 deg from the upper to the lower row 
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Fig. 10. Simulations results for a^—20° and fip — 30. We have selected different regions in the face on views to 
identify the locus of the different structures in the Zii-diagram (using the same color code in the face on view and the 
/•y-diagram) . Solid and dashed lines, circles and triangles are the same that in the previous figures. The dot-dashed 



lines in the lower ^w-diagram are the fits to the near and far 3-kpc arms as given by Dame & Thaddeus ( 2008 ) 
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Fig. 11. Simulations results for = 20°, % = 30 and q;^=90°, 75°, 60°, and 45° (from top to bottom). 
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Fig. 13. Simulations results for = 20°, fip — 30 and — 75° (with nuclear bar mass of 25% of bulge mass) 
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FlCE Ctl 




Fig. 14. Simulations results for = 20° , Hp = 30 and = 75° . The upper panel corresponds to the best fit to 
the 2MASS data with a lopsided bar (xi = 23 pc). The middle panel is a simulation with the same nuclear bar but 
Xi — 90 pc and the lower panel is a simulation with the same nuclear bar but Xi — —9 
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Table 3. Results of the fits to the star counts data of Fig. [5| using disc, bulge and nuclear bar. The following parameters 
have been taken from Tableland kept fix: = 10, = 1.117, = 0.513, = 0.362 K = 2.5, = 0.153. The 
bar angle is also fixed, the other parameters are free and used to obtain the best fit to the data. The first set of 
results have been obtained with xi fixed to 0. In the second set of results xi is a free parameter to fit. 
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Table 4. Same that Table [s] but with = 20, 0.884, yj^ = 0.489, ^ 0.375 K 2.5, h, 0.156. The first 

set of results have been obtained with xi fixed to 0. In the second set of results xi is a free parameter to fit. 
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Table 5. Same that Table [s] but with = 35, xE = 0.774, yE = 0.429, = 0.382 hr = 2.5, = 0.158. The first 
set of results have been obtained with xi fixed to 0. In the second set of results xi is a free parameter to fit. 
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350 


0.608E- 


02 





252E- 


03 





188 





126 


0.095 





141 


135.0 


0.030 





348 


0.609E- 


02 





252E- 


03 





170 





126 


0.094 





157 


120.0 


0.025 





348 


0.610E- 


02 





252E- 


03 





162 





122 


0.093 





171 


105.0 


0.022 





348 


0.610E- 


02 





252E- 


03 





163 





119 


0.091 





178 


90.0 


0.022 





348 


0.611E- 


02 





252E- 


03 





164 





124 


0.089 





171 


75.0 


0.022 





349 


0.611E- 


02 





252E- 


03 





173 





127 


0.091 





158 


60.0 


0.025 





349 


0.609E- 


02 





252E- 


03 





193 





126 


0.095 





140 


45.0 


0.031 





348 


0.609E- 


02 





252E- 


03 





222 





137 


0.095 





112 


30.0 


0.044 





348 


0.606E- 


02 





252E- 


03 





308 





149 


0.100 





073 


15.0 


0.089 





351 
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Table 6. Same that Table [s] but with = 45, 0.693, yjf = 0.367, = 0.385 K ^ 2.5, h, ^ 0.159. The first 

set of results have been obtained with xi fixed to 0. In the second set of results xi is a free parameter to fit. 



Bulge 




Disk 












Bar 
























h 






6 














2 


pE 




pE 




X 





Vo 


^0 


Po 


a*" 




X 














kpc 


kpc 


kpc 






deg 


kpc 






0.761E- 


02 





255E- 


03 





211 





167 


0.105 





095 


165.0 











389 


0.763E- 


02 





255E- 


03 





204 





154 


0.103 





106 


150.0 











382 


0.766E- 


02 





255E- 


03 





184 





143 


0.099 





126 


135.0 











378 


0.767E- 


02 





256E- 


03 





173 





138 


0.097 





141 


120.0 











376 


0.767E- 


02 





255E- 


03 





170 





136 


0.097 





147 


105.0 











375 


0.767E- 


02 





255E- 


03 





168 





136 


0.096 





149 


90.0 











376 


0.767E- 


02 





255E- 


03 





170 





136 


0.097 





147 


75.0 











376 


0.767E- 


02 





256E- 


03 





172 





138 


0.097 





141 


60.0 











377 


0.766E- 


02 





255E- 


03 





183 





143 


0.099 





126 


45.0 











379 


0.763E- 


02 





256E- 


03 





204 





151 


A 1 AO 

U.1U2 





107 


30.0 











384 


0.761E- 


02 





256E- 


03 





217 





161 


0.106 





092 


15.0 











391 


0.764E- 


02 





257E- 


03 





204 





153 


0.102 





099 


0.0 











390 


0.766E- 


02 





255E- 


03 





288 





150 


0.098 





076 


165.0 





079 





350 


0.767E- 


02 





255E- 


03 





220 





134 


0.096 





112 


150.0 





042 





347 


0.768E- 


02 





255E- 


03 





186 





126 


0.096 





141 


135.0 





030 





347 


0.770E- 


02 





255E- 


03 





166 





134 


0.092 





151 


120.0 





025 





345 


0.771E- 


02 





255E- 


03 





163 





122 


0.090 





171 


105.0 





022 





346 


0.771E- 


02 





255E- 


03 





160 





119 


0.090 





178 


90.0 





022 





346 


0.771E- 


02 





255E- 


03 





165 





122 


0.090 





170 


75.0 





022 





346 


0.770E- 


02 





255E- 


03 





171 





122 


0.092 





162 


60.0 





025 





346 


0.769E- 


02 





255E- 


03 





193 





126 


0.094 





140 


45.0 





031 





346 


0.768E- 


02 





255E- 


03 





227 





133 


0.095 





112 


30.0 





045 





347 


0.765E- 


02 





256E- 


03 





308 





146 


0.100 





073 


15.0 





090 





350 



